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Abstract 

In this paper we show existence and uniqueness of a solution for a system of m variational partial 
differential inequalities with inter-connected obstacles. This system is the deterministic version of 
the Verification Theorem of the Markovian optimal m-states switching problem. The switching cost 
functions are arbitrary. This problem is in relation with the valuation of firms in a financial market. 
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1 Introduction 

In this work we are concerned with the following system of m variational partial differential inequalities 
with inter-connected obstacles: 

m.in{vi(t,x) — max (— gij (t, x) +Vj(t,x)), 

-d t Vi(t, x) - Avi(t, x) - ipi(t, x)} = 0, V(t, x) G [0, T] x M k , i£l= {1, m}, 
Vi(T,x) = 0. 

where g^, ipi are continuous functions, A an infinitesimal generator associated with a diffusion process 
and finally I~ l := {1, i — 1, % + 1, m}. 

This system is the deterministic version of the Verification Theorem of the optimal multi-modes 
switching problem in finite horizon. This problem, of real option type, can be introduced with the help 
of the following example: 
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Assume we have a power station/plant which produces electricity and which has several modes of 
production, e.g., the lower, the middle and the intensive modes. The price of electricity in the market, 
given by an adapted stochastic process (Xt)t<T, fluctuates in reaction to many factors such as demand 
level, weather conditions, unexpected outages and so on. On the other hand, electricity is non-storable, 
once produced it should be almost immediately consumed. Therefore, as a consequence, the station 
produces electricity in its instantaneous most profitable mode known that when the plant is in mode 
i Gl, the yield per unit time is given by means of ipi and, on the other hand, switching the plant from 
the mode i to the mode j is not free and generates expenditures given by gtj and possibly by other 
factors in the energy market. So the manager of the power plant faces two main issues: 

(i) when should her decide to switch the production from its current mode to another one? 

(ii) to which mode the production has to be switched when the decision of switching is made? 

In other words she faces the issue of finding the optimal strategy of management of the plant. This 
issue is in relation with the price of the power plant in the energy market. 

For decades, optimal switching problems have attracted a lot of research activity (see e.g. [H [51 
O El El El HOl HD H21 [13 [13 [HI HH EDI [22I E3l ED ESI ESI [HI ISH [29] and the references therein). 
Especially in connection with valuation of firms and, questions related to the structural profitability 
of investment project or an industry whose production depends on the fluctuating market price of a 
number of underlying commodities or assets, ... Several variants of the problem we deal with here 
have been considered. In order to tackle those problems, authors use mainly two approaches. Either 
a probabilistic one ( [101 [HI [HI CEH1 E21 [23] ) or an approach which uses partial differential inequalities 
(PDIs for short) ([B [2j H3J W\ US EE]). 

The PDIs approach turns out to study and to solve, in some sense, the system of m PDIs with 
inter-connected obstacles (jl.ip . Amongst the papers which consider the same problem as ours, and 
in the framework of viscosity solutions approach, the most elaborated works are certainly the ones by 
Tang and Yong [26J, on the one hand, and by Djehiche et al. [11 j . on the other hand. In [26J, the 
authors show existence and uniqueness of a solution for (|l.ip . Nevertheless the paper suffers from two 
facts: (i) the growth exponent at infinity of the functions ipi should be smaller that 2 ; (ii) the switching 
cost functions gij should not depend on x. The first issue of has been treated by Djehiche et al.|ll| 
since in their paper the authors show existence of the solution for 11.11 in the case when the growth of the 
functions tpi is of arbitrary polynomial type. The second issue of [26] , i.e. considering the case when 
depending also on x, was right now, according to our knowledge, an open problem. Note that in [llj . 
the question of uniqueness is not addressed. Therefore the main objective of our work, and this is the 
novelty of the paper, is to show existence and uniqueness of a solution in viscosity sense for the system 
when the functions ipi and g^ are continuous depending also on x and satisfy an arbitrary polynomial 
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growth condition. We show also that the solution is unique in the class of continuous functions with 
polynomial growth. 

This paper is organized as follows: 

In Section 2, we formulate the problem and we give the related definitions. In Section 3, we introduce 
the optimal switching problem under consideration and give its probabilistic Verification Theorem. It 
is expressed by means of a Snell envelope of processes. Then we introduce the approximating scheme 
which enables to construct a solution for the Verification Theorem. Moreover we give some properties 
of that solution, especially the dynamic programming principle. Section 4 is devoted to the connection 
between the optimal switching problem, the Verification Theorem and the system of PDIs (jl.ip . This 
connection is made through backward stochastic differential equations with one reflecting obstacle in 
the case when randomness comes from a solution of a standard stochastic differential equation. Further 
we provide some estimate for the optimal strategy of the switching problem which in combination 
with the dynamic programming principle plays a crucial role. Finally we show that system (11. ip has a 
solution. In Section 5, we show that the solution of (jl.ip is unique in the class of continuous functions 
which satisfy a polynomial growth condition. □ 

2 Assumptions and formulation of the problem 

Throughout this paper T (resp. k) is a fixed real (resp. integer) positive constant. Let us now consider 
the followings: 

(i) b : [0, T] x M k -» M k and a : [0, T] x M k -> M kxd are two continuous functions for which there 
exists a constant C > such that for any t G [0, T] and x, x' G M k 

\b(t,x)\ + \<x(t,x)\ < C(l + \x\) and \a(t, x) - a(t, x')\ + \b(t, x) - b(t, x')\ < C\x - x'\ (2.1) 

(ii) for i,j £l = {1, m}, gij : [0, T] x M k — > M and ifii : [0, T] x ]R k — > 1R are continuous functions 
and of polynomial growth, i.e., there exist some positive constants C and 7 such that for each i, j G I: 

\Mt,x)\ + \9ij{t,x)\<C(l + \x\i), V(t,x) G [0,T}xM k . (2.2) 

Moreover we assume that there exists a constant a > such that for any (t,x) G [0, T] x M k , 

min{gij(t,x),i,j G X, i ^ j} > a. (2.3) 

We now consider the following system of m variational inequalities with inter-connected obstacles: 
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V« el 

min <^ Vi(t,x) - max (-gij(t,x) + Vj(t,x)), -d t Vi(t,x) - Avi(t,x) - ipi(t,x) > = 0, 

I iex- 1 J (2.4) 

Vi(T,x) = 0, 

where I~ l := I — {i} and ^4 is the following infinitesimal generator: 

a- \ e (-*)«('-)s^ + £ ; (2 - 5) 

i,j=l,k J i=l, fc 

hereafter the superscript (*) stands for the transpose, TV is the trace operator and finally < x,y > is 
the inner product of x, y G iR fc . 

The main objective of this paper is to focus on the existence and uniqueness of the solution in 
viscosity sense of (|2.4p whose definition is: 



Definition 1 Let (v±, ...,v m ) be a m-uplet of continuous functions defined on [0, T] x M k , M-valued 
and such that such that Vi(T,x) = for any x G M k and i £l. The m-uplet (v\, ...,v m ) is called: 



(i) a viscosity super solution (resp. subsolution) of the system \2.J$ if for each fixed i G X, for any 
(toj^o) ^ [0>^1 x lR k an d an y function ipi G C 1,2 ([0,T] x M k ) such that ipi(to,xo) = Vi(to,xo) and 
(to,xo) is a local maximum of ipi — v i (resp. minimum), we have: 



mini Vi(t ,x ) - max (-gij(t ,x ) +Vj(t ,x )), 

I jex~ i (2.6) 

-d t (pi(to,x ) -A<pi(t ,x ) -ipi(t ,x )} > (resp. < 0). 
(ii) a viscosity solution if it is both a viscosity supersolution and subsolution. □ 

There is an equivalent formulation of this definition (see e.g. [6]) which we give because it will be 
useful later. So firstly we define the notions of superjet and subjet of a continuous function v. 

Definition 2 Let v G C((0, T) x M k ), (t,x) an element of (0, T) x M k and finally Sk the set of k x k 
symmetric matrices. We denote by J 2 ' + v(t, x) (resp. J 2, ~v(t,x)), the superjets (resp. the subjets) of v 
at (t,x), the set of triples (p,q,X) G 1R x ]R k x such that: 

v{s, y) < v{t, x) + p(s -t) + (q,y - x) + \(X(y -x),y-x) + o(\s - t\ + \y - x\ 2 ) 
(resp. v(s, y) > v(t, x) + p(s — t) + (q,y — x) + \(X(y — x), y — x) + o(\s — t\ + \y — x\ 2 )).U 

Note that if ip — v has a local maximum (resp. minimum) at (t, x), then we obviously have: 



(D t (p(t, x),D x tp(t, x), Dl x <p(t, x)) G J 2 ' v(t,x) (resp. J 2 ' + v(t, x)). 



□ 



We now give an equivalent definition of a viscosity solution of the parabolic system with inter- 
connected obstacles (12.41). 



Definition 3 Let (v\, v m ) be a m-uplet of continuous functions defined on [0, T] x M k , M-valued 
and such that (v±, v m )(T, x) = for any x G M k . The m-uplet (v±, ...,v m ) is called a viscosity 
supersolution (resp. subsolution) of \2.J$ if for any i £l, (t,x) G (0,T)xM k and (p, q, X) G J 2, ~Vi(t,x) 
(resp. J 2,+ Vi(t,x)), 

min l Vi(t, x) — max i—gijit, x) + Vj(t, x)), —p — — Trier* Xa] — (b, q) — ibAt, x) 1 > (resp. < 0). 
I i€X-« 2 J 

It is called a viscosity solution it is both a viscosity subsolution and supersolution .□ 

As pointed out previously we will show that system (|2,4p has a unique solution in viscosity sense. 
This system is the deterministic version of the optimal m-states switching problem which is well docu- 
mented in [11] and which we will describe briefly in the next section. 

3 The optimal m-states switching problem 
3.1 Setting of the problem 

Let (£l,J-,P) be a fixed probability space on which is defined a standard d-dimensional Brownian 
motion B = (B t )o<t<T whose natural filtration is {J^ := a{B s ,s < t})o<t<T- Let F = (F^o^t^T be 
the completed filtration of (J r t °)o<t<T with the P-null sets of hence (J-t)o<t<T satisfies the usual 
conditions, i.e., it is right continuous and complete. Furthermore, let: 

- V be the cr-algebra on [0, T] X f2 of F-progressively measurable sets; 

- A4 2 ' k be the set of "P-measurable and iR fc -valued processes w = (wt)t<T such that E[Jq \w s \ 2 ds] < 
oo and S 2 be the set of P-measurable, continuous processes w = [w t )t<T such that £'[sup t<T \ w t \ 2 \ < oo; 

- for any stopping time r G [0, T], T T denotes the set of all stopping times 9 such that r < 9 < T. 

□ 

The problem of multiple switching can be described through an example as follows. Assume we 
have a plant which produces a commodity, e.g. a power station which produces electricity. Let I be the 
set of all possible activity modes of the production of the commodity. A management strategy of the 
plant consists, on the one hand, of the choice of a sequence of nondecreasing stopping times (r n ) n >i 
(i.e.T n < T n+ \ and to = 0) where the manager decides to switch the activity from its current mode 
to another one. On the other hand, it consists of the choice of the mode £ n , a r.v. T Tn -measurable 
with values in X, to which the production is switched at r n from its current mode. Therefore the 
admissible management strategies of the plant are the pairs (5, £) := ((x„) n >i, (^ n ) n >i) and the set of 
these strategies is denoted by V. 

Let now X := (A"t)o<t<T be an P-measurable, iR fc -valued continuous stochastic process which stands 
for the market price of k factors which determine the market price of the commodity. On the other 
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hand, assuming that the production activity is in mode 1 at the initial time t = 0, let {ut)t<T denote 
the indicator of the production activity's mode at time t G [0, T] : 

U t = l [0tTl] (t) + J2^ {Tn ,r n+1 ](t). (3.1) 
n>l 

Then for any t <T, the state of the whole economic system related to the project at time t is represented 
by the vector : 

(t,X t ,u t ) G [0,T] x M k xl (3.2) 

Finally, let ipi(t,Xt) be the instantaneous profit when the system is in state (t,Xt,i), and for 
i 7^ j, let gij(t, Xt) denote the switching cost of the production at time t from current mode % 
to another mode j. Then if the plant is run under the strategy (5, £) = ((x„)„>i, (£ n ) n >i) the expected 
total profit is given by: 

J(S,0 = E[ / ip Us (s,X s )ds - YV UTn (r n ,X Tn )l [Tn<T] ]. 

J n >l 

Therefore the problem we are interested in is to find an optimal strategy, i.e, a strategy (<5*,£*) such 
that J{5*,C) > J(S,0 for any (6,0 G V. 

Note that in order that the quantity J(5, £) makes sense we assume throughout this paper that for 
any i,j G X the processes (ipi(t, X t )) t <T and (gij(t, X t )) t <T belong to M 2 ' 1 and S 2 respectively. On the 
other hand there is a bijective correspondence between the pairs (<5, £) and the pairs (<5, u). Therefore 
throughout this paper one refers indifferently to (5, £) or (8, u). 

3.2 The Verification Theorem 

To tackle the problem described above Djehiche et al. [11] have introduced a Verification Theorem 
which is expressed by means of Snell envelope of processes. The Snell envelope of a stochastic process 
{i]t)t<T of S 2 (with a possible positive jump at T) is the lowest supermartingale R(rj) := (R(r])t)t<T of 
S 2 such that for any t <T, R(T])t > i]t- It has the following expression: 

Vt < T,R{rj) t = esssup Te T t E[rj T \F t ] and satisfies R{j])t = rjr- 

For more details on the Snell envelope notion on can see e.g. [4, ll"l I16j. 

The Verification Theorem for the m-states optimal switching problem is the following: 

Theorem 1 ( 11 Jj /. Th.l) Assume that there exist m processes (Y i := (Y^)o<t<T)2 = l,...,m) of S 2 
such that: 

Vt < T, Y* = ess sup T > t E[f t T ^(s, X s )ds + max (-^-(r, X T ) + F/)1 [t<t] \T t ], Y* = 0. (3.3) 

j"6Z-* 

T/ien: 
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(i) Yq = sup J(S,u). 

(<5,C)ex> 

(ii) Define the sequence ofF-stopping times 5* = (r^) n >i as follows : 
t{ = inf{s > 0, Y* = max (-<?y (s, X s ) + Y/)} A T, 

t* = inf^r*^, max (-<?„ fc (a, X s ) + Y s fc )} A T, for n > 2, 

fcgx\{«_. } ""I 

n — 1 



where: 



«rf ~ E jl { max (- Jlfc (T 1 *,X .)+^)=- mj (Tf,X .)+y;,}' 
jG-t fcei- 1 1 i 1 i 

/or any n > 1 and i > r*, Y t " T ™ = X) 1[« *=i]Y/ 
for any n > 2, u T * = / on the set 



I t ,r^ aX ^ ( ~^* A<> X <)+ Y k) = ~9u T , (r*,X T *) + Yl, \ 

\k£l\{u T * } n-1 n n-1 "J 

with g u k (T*,X T *) = E 1[ =i]ffifc«,X r *) andX\{u r * } = £ ![«. =j] x-i - 

n-1 j-gj «-l ' jgj n-1 

T/ien the strategy (5*,u*) is optimal i.e. J(5*,u*) > J(5,u) for any (5,u) £ T>. □ 

The issue of existence of the processes Y , ...,Y m which satisfy (|3.3[) is also addressed in [11]. Also 
for n > let us define the processes (Y" ,n , Y m ' n ) recursively as follows: for i € I we set, 

Y/>° = ess sup T > t E[£ iJji(s,X s )ds\Ft], < t < T, (3.4) 

and for n > 1, 

y/' n = ess sup T>t E[ [ T A(s,X s )ds+ m a x(-g lk ( T ,X T ) +Y T k > n - 1 )l [r<T] \f t ], < t < T. (3.5) 

Jt k<=l- 1 

Then the sequence of processes ((Y ' n , Y m ' n )) n >o have the following properties: 

Proposition 1 (111], Pro. 3 and Th.2) 

(i) for any i G X and n>0, the processes Y 1,n , ...,Y m ' n are well-posed, continuous and belong to S 2 , 
and verify 

rT 

Vt < T, y/' n < y/' n+1 < E[ / {max {^(s, X s )\}ds\F t ]; (3.6) 

Jt l=1 > m 

(ii) there exist m processes Y 1 , Y m of S 2 such that for any i £ I: 
(a) Vt < T, y/ = limn-.oo / Y t l ' n and 

E[sup \Y l s ' n -y?| 2 ] -> as n -> +oo 



(6) Vt < T, 

Yl = ess sup T>t E[ [ iPi(s, X s )ds + max (-g ifc (r, X T ) + ^Wt] l^t] (3-7) 

i.e. Y" , ...,Y m satisfy the Verification Theorem^; 
(c) Vt < T, 

f T 

E[ ^ Us (s,X s )ds-Y^ ^r„)l[r„<T|]|^] (3-8) 

Jt n>l 

where T>\ = {(<5, £) = ((r n ) n >i, (£n)n>i) such that uq = i and t\ > t}. This characterization 
means that if at time t the production activity is in its regime i then the optimal expected 
profit is Y t l . 

(of) the processes Y , Y m verify the dynamical programming principle of the m-states optimal 
switching problem, i.e., Vi < T , 

iP U3 (s,X s )ds- 9u Tk _ l u Tk {T k ,X Tk )\ [Tk<T] + \ [Tn<T] Y^\T t ].n 

l<k<n 

(3-9) 

Note that except (ii — d), the proofs of the other points are given in [11] . The proof of (ii — d) can 
be easily deduced in using relation (|3.7p . Actually from (|3.7j) for any i El, t E [0,T] and (5, £) E T>\ 
we have: 

Y}>E[j T \ Ua {s,X s )ds- Y, K-iH^' 1 ^ 1 ^^!^^" 1 ^- (3 ' 10) 

Jt l<k<n 

Next using the optimal strategy we obtain the equality instead of inequality in f|3. 10]) . Therefore the 
relation (|3.9|) holds true. □ 

Remark 1 Note that the characterization \3. 8]) implies that the processes Y , Y m of S 2 which satisfy 
the Verification Theorem are unique. 

4 Existence of a solution for the system of variational inequalities 

4.1 Connection with BSDEs with one reflecting barrier 

Let (t, x) E [0, T] x M k and let {X l s x ) S <T be the solution of the following standard SDE: 

dXf = b(s, Xl x )ds + a(s, Xf )dB s for t < s < T and Xf = x for s < t (4.1) 

where the functions b and a are the ones of ()2. 1 j) . These properties of a and b imply in particular that 
the process (Xl x )o< s <T solution of the standard SDE (|4.ip exists and is unique, for any t € [0, T] and 
x E M k . 



The operator A that is appearing in (|2.5p is the infinitesimal generator associated with X t,x . In the 
following result we collect some properties of X t,x . 

Proposition 2 (see e.g. \2J$ ) The process X tx satisfies the following estimates: 
(i) For any q > 2, there exists a constant C such that 

E[ sup \X tx \ q ] < C(l + \x\ q ). (4.2) 

0<s<T 

(ii) There exists a constant C such that for any t,t' G [0, T] and x,x' G M k , 

E[ sup \X tx - X* V | 2 ] < C{1 + |x| 2 )(|x - x'\ 2 + \t- t'\).a (4.3) 

0<s<T 

We are going now to introduce the notion of a BSDE with one reflecting barrier introduced in |15j . 
This notion will allow us to make the connection between the variational inequalities system (|2.4D and 
the m-states optimal switching problem described in the previous section. 

So let us introduce the deterministic functions / : [0, T] x ]R k+1+d —> M, h : [0, T] x M k —> M and 
g : M k — > ]R continuous, of polynomial growth and such that h(x, T) < g(x). Moreover we assume that 
for any (t,x) £ [0, T] x M k , the mapping (y,z) € M 1+d i— > f(t,x,y,z) is uniformly Lipschitz. Then we 
have the following result related to BSDEs with one reflecting barrier: 

Theorem 2 (115$ , Th. 5.2 and 8.5) For any (t,x) E [0, T] x M k , there exits a unique triple of processes 
(Y t ' x ,Z t ^,K t ' as ) such that: 

Y tx , K tx G S 2 and Z tx G M 2 ' d ; K tx is non- decreasing and Kjf = 0, 
< Y s tx = g{X%) + Jf f(r, X tx X x ^ Z tx )dr - jf Z tx dB r + K$ - K tx , s<T (4.4) 
k Y tx > h(s,X tx ), Vs < T and f^(Y r tx - h(r,X tx ))dK tx = 0. 

Moreover the following characterization of Y t,x as a Snell envelope holds true: 

Vs < T, Y*>* = esssup Te r t E[f t f(r, X* x ,Y r tX : + h(r, X tx )\ [r<T] + g(X^)l [r=T] (4.5) 

On the other hand there exists a deterministic continuous with polynomial growth function u : 
[0, T] x M k -> M such that: 

y s e[t,T],Y s t>x = <s,X^). 

Moreover the function u is the unique viscosity solution in the class of continuous function with poly- 
nomial growth of the following PDE with obstacle: 

min{u(t, x) — h(t, x), —dtu(t, x) — Au(t, x) — f(t, x, u(t, x), a(t, x)*Vu(t, x))} = 0, 
u(T, x) = g{x).n 
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4.2 Existence of a solution for the system of variational inequalities 

Let (y/- M ,...,y, m ' t3 )o<s<T be the processes which satisfy the Verification Theorem [1] in the case when 
the process X = X l ' x . Therefore using the characterization (|4.5p . there exist processes K l ' tx and Z l ' tx , 
i El, such that the triples (Y htx , Z l ' tx , K l ' tx ) are unique solutions (thanks to Remark[T|) of the following 
reflected BSDEs: for any i = 1, ...,m we have, 

' Y i,tx^ K i,tx e S 2 and Z i,tx e M 2,d. K i,tx is non .decreasing and K l ' tx = 0, 

Y l s ' tx = J S T Vi(r, J^)du - Z l r ' tx dB r + - K^ tx , < s <T, Y^ tx = 0, 
< te > max (- 9ij (s, Xt x ) + Yj> tx ), < a < T, ^ 

J T (y;- te - max(-^(r,X^) + Y r j ' tx ))dKp tx = 0. 
Moreover we have the following result. 

Proposition 3 There are deterministic functions v 1 , ...,v m : [0, T] x ]R k — » iR suc/i i/mi: 

V(t,x) G [0,T] x iR fe ,Vs G [t,T],y?' te = v\s,X? ), i = l,...,m. 
Moreover the functions v l , i = 1, ...,m, are /ouier semi- continuous and of polynomial growth. 

Proof: For n > let (y"' 1 ' , Y s n ' m ' tc °)o<s<T be the processes constructed in f|3.4[) - (|3. 5j) . Therefore 
using an induction argument and Theorem 2 there exist deterministic continuous with polynomial 
growth functions v^ n (i = l,...,m) such that for any (t,x) G [0,T] x M k , Vs G [t,T], Y^' tx = 
v l ' n (s, X tx ). Using now inequality (|3.6p we get: 

y n,i,tx < ^H-l,i,te < /" T { max |^( S) 

Ji »=l,m 

since Y™ ,l,tx is deterministic. Therefore combining the polynomial growth of ipi and estimate (|4.2p for 
X'^ we obtain: 

?/' n (t,x) < v^ n+1 {t,x) < C(l + |x| p ) 

for some constants C and p independent of n. In order to complete the proof it is enough now to set 
:= lim^oo a?), G [0,T] x since Y^ tx / Y^ tx as n -» cx). □ 

We are now going to focus on the continuity of the functions u 1 , ...,v m . But first let us deal with 
some properties of the optimal strategy which exist thanks to Theorem 1. 

Proposition 4 Let (S,u) = ((T n ) n >i, (£ n ) n >i) be an optimal strategy, then there exist two positive 
constant C and p which does not depend on t and x such that: 

Vn > 1, P[r n < T] < C{l + \ X \ P \ (4. 7 ) 

n 
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Proof: Recall the characterization of (|3.8p that reads as: 

Yq M = sup {s ,u)evE[fQ ip Ur (r, X^dr - ^2 k >i 3^ k _ x u Tk (r k , X% )% k< T}\- 
Now if (5,u) = ((r n ) n >i, (£ n )n>i) is the optimal strategy then we have: 

Y^ tX = E[f^ Ur ( r ,X?)dr - Ek>i 9u Tk _ lUrk (n, X^ k )l [Tk<T] }. 
Taking into account that gij > a > for any i / j we obtain: 

E[Zk=i,n «l[r fc <T]] + Y^ tx < E[fo ^u r (r, X*)dr - £ fc > n+1 9u Tk _ x u Th (r fc , X^)l [rfc<T] ]. 
But for any k < n, [r n < T] C [r^ < T] then: 

anP[r n < T]+Y^ tx < E[ff ^(r, X tx )dr - £ fc > n+1 g UTk _,u Tk (r k , X tx )\ [Tk<T] ] 



<E[J^ Ur (r,X^)dr]. 



and then 



na 



P[r n < T] < E[J* | W (r, X?) \ dr] - Y c 



l,tx 



<E[ff\^ Ur (r,X^)\dr]-Y^ tx . 
Finally taking into account the facts that tpi and Y l,0,tx are of polynomial growth and estimate (|4.2p 
for X tx to obtain the desired result. Note that the polynomial growth of Y 1,0,tx stems from Proposition 
3. □ 



Remark 2 The estimate \4. 7\ ) is also valid for the optimal strategy if at the initial time the state of 
plant is an arbitrary i G 1. □ 



Next for i G X let (yl' tx , zl' tx , kl' tx ) s <T be the processes defined as follows: 



yi,tx^i,tx e £2 an( j z i,tx g _A/fAd. /.i,ta: j g non-decreasing and A: 
i,tx I „/. /_ v-ta;\ii . . j„ / -I-i-.it> . i. iJ.'->- 



i,tx 



^i(r,X™)l [r > t] dr - / zJ? te dB r + - kf x ', < s < T, $ 



yjte > := max{ _ 3 ( t v 8 ,X$° a ) + yi' tx )}, Vs < T, 
J > tx lf x )dk\: tx = 0. 



(4.8) 



K JO 



The existence of (y l ' tx , z htx , k t,tx ), i <G X, is obtained in the same way as the one of (y htx , Z l ' tx , K t,tx ). 
On the other hand, thanks to uniqueness (see once more Remark [I]), for any (t,x) E [0,T] x M k , for 



any s G [0, t] we have yl' tx = Y t 



A,tx i,tx 
t ' z s 



and k 



i,tx 



0. 



We are now ready to give the main Theorem of this article. 

Theorem 3 The functions (v 1 , ...,v m ) : [0,T]xM k -> 1R are continuous and solution in viscosity sense 
of the system of variational inequalities with inter- connected obstacles l[2.4\)- 
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Proof: First let us focus on continuity and let us show that v 1 is continuous. The same proof will be 
valid for the continuity of the other functions v l (i = 2, ...,m). 
First the characterization (|3.8p implies that: 



yl' tx = sup E[ [ 4> Us (s, Xl x )l [s > t] ds - V g UTn lUrn (r n V t, X% vt )l [Tn<T] } 



On the other hand an optimal strategy (5* , ) exists and satisfies the estimates (I4.7P with the same 
constants C and p. Next let e > and (t',x') E B((t,x),e) and let us consider the following set of 
strategies: 

C(\ + (f + \x\) p ) 

D := {(6,0 = ((r n )„>i, (£„.)„>o) G ^ such that Vn > l,P[r n < T] < 1 V 1 M ; }. 

~~ ~~ n 

Therefore the strategy (<5*,£*) belongs to D and then we have: 

vl' tx = sup me6 E[fQ^ Us (s,Xt x )TL[ s >t]ds - E n >i ^ V ^v^Ir^T]] 
= s u P(s,u)eD E Uo n ^u s {s,Xt x )l [s > t] ds 

~ Ei<fe< n 9u Tk _ x u Tk (t y T k ,X^ Tk )l [Tk<T] + l[ r , l<T ]yr,r n ' te ] 

and 

2/0'* V = SUp m( - D E[f* lp Us (S, X* V )l[s>i']<fe - E n >l ft*r n _i«^, ( T « V ^rlvt')l[r„<T]] 

= sup^g^Pf/o 7 "" ^ Ua (s,X^ x ')\ {s > tl] ds 

The second equalities follow from the dynamical programming principle. It follows that: 

y i,tv _ y i,t x < ^^^[/^{^(s.X^)!^] - i> U3 (s,Xl x )l [s > t] }ds 

(4.9) 

+ ^[r n <T]{yr T y X ' -Vr T n n ' tX }] 

Next w.l.o.g we assume that t' < t. Then from (14. 9j) we deduce that: 

- EkkJ^h^J*' V T k,X^ Tk ) - 9u Tk ^u Tk (t V T k ,X$ Tk )}l [Tk<T] 



^[r n <T]{Vr T y X -yZ n,tX }] 



^ E [Jo max j=l,m\i>j(s,Xt' x ') -i/}j(s,Xl x ))\l [s > t] +max j=ltm \il> j (s,Xl' x ')\l [t/ < :S<t] }ds 
+n max^ jgJ {sup s < T \gij(t' V s, X$ x ' s ) - g^it V s, X$ s )\} 
+^P(6,n)eD(P[rn < T\)t {2E[{y u T ;r t,x ') 2 + (y^' te ) 2 ])i 



(4.10) 

In the right-hand side of (|4.10p the first term converges to as (t',x') — > (t,x). Next let us show that 
for any ij^jEl, 

E[sup \ 9ij (t' V s, X^' s ) - 9ij (t V s, Xj x s )\] ^ as [t' , x') -» (t, x). 

s<T 
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Actually for any g>0we have: 

\ gij (t' -^(^8,^)1 < \ gij tfV8,X$*)- gij (t V8,X$*)\1 r|JftV ,<„!+ 

|^(f V s, - 0y(t V s, X*'^)|l n w,/ (> , + \ 9ij (t V a, X$*' s ) - gij (t V s, X$ a )\. 

Therefore we have: 

£[sup s < T \ 9lJ (t' V s, Xf-'J - gij {t V s, X$ s )\] < 

E[sup s < T {\ gij {t' V s, X$*j a ) - gij (t V s, X$Q\l llxt y ,< ,}]+ 
£[sup s < T {| fti (i' V s, X$*) - 9ij (t V a, X™)\}l [mpsiT |> c] ]+ 
E[sup s < T {\ 9ij (t V s, X$* a ) - 9ij (t V s, X|^)|}l [sup ^ T |x^'j +S u Ps < T m> B ]]+ 
E[sn Vs < T {\ gij (t V a, Xffi s ) - 9lj (t V s, X$ 8 )\}l [Bap ^ T \ x t>*>j +sup ^ T \ x t* s] < g] ] 

But since 9 ij is continuous then it is uniformly continuous on [0, T] x {x € iR fc , |x| < f?}. Henceforth for 
any ei > there exists rj ei > such that for any \t — t'\ < r] ei we have: 



su V {\ 9ij (t' Vs,!^) - W (t Vs,X^J|l [|xtV|< ,} < ei. (4.11) 

Next using Cauchy-Schwarz's inequality and then Markov's one with the second term we obtain: 

E[ S up{\ gij (t'Vs,Xt:*' s ) - 9ij (t V8,xtf a )\}l. <T ix^'\>J < C (! + W\ P )e~* (4-12) 

where C and p are real constants which are bound to the polynomial growth of 9 ij and estimate (|4.2p . 
In the same way we have: 

£[sup{| 5ii (iV S ,X*;*:)-^ (4.13) 



Finally using the uniform continuity of cjij on compact subsets, the continuity property (14.3(1 and the 
Lebesgue dominated convergence theorem to obtain that 

E[sM\ 9ij {t V a, X^ s ) - 9ij (t V s, A^J|lL [sups<T \x?*' 3 \+ S u Ps < T \x&.\< B ]i ~> as x ') (*' x )- ( 4 - 14 ) 
Taking now into account (|4.1ip - (|4.14p we have: 

limsup £[sup \ gij (t' V s, X$* 8 ) - 9ij {t V s, X^ s )\] < ei + C(l + 

ft' ,2') S < T 

As ei and £ are arbitrary then making e\ — * and g — > +oo to obtain that: 

lim £[sup \ 9ij (t'Vs, X^ s ) - 9ij (t V s, X%,)\] = 0. 
(1/ ,x')-*{t,x) S <T 

Thus the claim is proved. 

Finally let us focus on the last term in (|4,10p . Since (8,u) G D then: 

™ Pis , u)e B(P[Tn < T])\{2E[{y u T r t,X 'f + (l^**) 2 ])* < sup (M65 (2i?[(^ V ) 2 + (y^) 2 ])^ 

< Cn"5(l + \ X \P+ \x'\p) 
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where C and p are appropriate constants which come from the polynomial growth of ipi , i € I, estimate 
(|4.2p for the process X tx and inequality (|3.6p . Going back now to (|4.10p . taking the limit as (t',x') — > 
(i, x) to obtain: 

limsup yg'* x < yo' tx + Cn~ 2 (1 + 2|xj p ). 

(t',a:')-»(t,a;) 

As n is arbitrary then putting n — > +00 to obtain: 

limsup y ' < y ' . 

(f ,x')->(t,x) 

It implies that: 

limsup yl' tx = Y^' tx = v l {t' ,x') < y^ tx = Y^' tx = x). 

Therefore v 1 is upper semi-continuous. But v 1 is also lower semi-continuous, therefore it is continuous. 
In the same way we can show that v 2 ,...,v m are continuous. As they are of polynomial growth then 
taking into account Theorem 2 to obtain that (v 1 ,...,t> m ) is a viscosity solution for the system of 
variational inequalities with inter-connected obstacles (|2.4p . □ 

5 Uniqueness of the solution of the system 

We are going now to address the question of uniqueness of the viscosity solution of the system (|2.4p . 
We have the following: 

Theorem 4 The solution in viscosity sense of the system of variational inequalities with inter- connected 
obstacles \2.J$ is unique in the space of continuous functions on [0, T] x R k which satisfy a polynomial 
growth condition, i.e., in the space 

C :={</?: [0, T] X M k — > M, continuous and for any 

(t,x), \ip(t,x)\ < C(l + |x| 7 ) for some constants C and 7}. 

Proof. We will show by contradiction that if u\, u m and W\, w m are a subsolution and a supersolu- 
tion respectively for (12. 4p then for any i = 1, m, Ui < Wi. Therefore if we have two solutions of (|2.4f) 
then they are obviously equal. Actually for some R > suppose there exists (t,x,i) G (0,T) x Br x Z 
(B R := {x G M k ; \x\ < R}) such that: 

uj(t,x) - wj(t,x) = n > 0. (5.1) 

Let us take 6, A and (3 € (0, 1] small enough, so that the following holds: 

' (3T< I 

29\x\ 2 ~<+ 2 < 3 

__ (5-2) 
-\w-(t,x) < g 

^ t ^ 5" 
14 



Here 7 is the growth exponent of the functions which w.l.o.g we assume integer and > 2. Then, for a 
small e > 0, let us define: 

$*(t, x, y) = Ui(t, x)-(l- X) Wi (t, y) - i|x - y| 27 - 9(\x\ w + |y | 2 ^+ 2 ) + ft - j. (5.3) 

By the growth assumption on Uj and u^, there exists a (to, ^o> 2/o> *o) G (0, T] x x Z? K x Z, such that: 

<^ (t ,£ ,yo) = max ^(^^y)- 

(t,x,y,t) 

On the other hand, from 2<j?*°(t , £oj 2/o) > 3>e°(*o> x Cb ^0) + ^e O (*0)2/0)2/o)) we nave 

2^ko - 2/o| 27 < (ui (t ,x ) - u io (t ,y )) + (1 - X)(w io (t ,x ) - w io (t ,yo)), (5.4) 

and consequently ^\xq — yo| 27 is bounded, and as e — > 0, |xo — yo\ — > 0. Since Ui and Wj are uniformly 
continuous on [0,T] x Br, then t^I^o — yo| 27 —> as e — > 0. 



Next let us show that to < T. Actually if to = T then, 

$l(lx,x)<<f>i°(T,x ,y ), 

and, 

u-(t,x) - (1 - AH(i,z) - 2#|xj 2 ^ +2 + ft - j < 0T - ±, 



since Ui (T, xq) = Wi (T,y ) = 0. Then thanks to (|5.ip we have, 

V < -Xio^t,x) + f3T + 26\x\ 2 ^ +2 + f 
7? < fr?. 

which yields a contradiction and we have to G (0,T). We now claim that: 



Indeed if 



«i (*o,a;o)- max {-g io j(to,x ) + Uj(to,x )} > 0. (5.5) 



«io(*o,^o) - max {-g ioj (to,xo) + Uj(t ,x )} < 



then there exists k £ 2 10 such that: 

«io(*0,»o) < -9i k(to,X ) +Uk(ta,X ). 

From the supersolution property of Wi (to,yo), we have 

u>io(*o,yo) > max (-gi j(to,y ) +Wj(t ,y )) 
j£l-*o 

then 

w io (to,yo) > -9iok(to,yo) + w k (t ,yo)- 
15 



It follows that: 



Ui (to,x ) - (1 - X)w io (t ,y ) - (u k {t ,x ) - (1 - X)w k (to,y )) < (1 - A)g iofe (to, yo) - 5iofc(*o, ^o)- 

Now since > a > 0, for every i 7^ j, and taking into account of (|5.3p to obtain: 

$*°(io,xo,yo) - $e(io,£o,yo) < -a\ + g io k(to,yo) - gi k(to,x ). 

But this contradicts the definition of io, since gi Q k is uniformly continuous on [0, T] x 5^ and the claim 
(1531) holds. 

Next let us denote 



¥> e (t, x,!/) = l|x - y| 27 + #(|x| 2 ^ 2 + |y| 27+2 ) - # + ~. 



(5.6) 



Then we have: 



D t cp e (t,x,y) = -(3- -p, 

D x <p e (t,x,y) = ^(x - y)\x - y\ 2 ^ 2 + 6{2 1 + 2)x\x\ 2 \ 
D v <p e (t, x, y) = -l{x - y)\x - y\ 2 ^ 2 + 9{2 1 + 2)y\y\ 2 \ 



B(t,x,y) = D 2 cp e (t,x,y) 



l l at(x,y) -ai(x,y)\ / a 2 (x) 



-ax(x,y) at(x,y) 



+ 



(5.7) 



a 2 (y), 



with ai(x, y) = ^\x - y\ 21 2 I + 7(27 — 2)(x — y)(x - y)*\x - y| 27 4 and 
^ 02 (x) = 0(2-/ + 2)\x\ 2 ^I + 26» 7 (2 7 + 2)xx*\x\ 2 ^- 2 . 

Taking into account (15. 5p then applying the result by Crandall et al. (Theorem 8.3, [6]) to the function 

u io (t,x) - (1 - X)w io (t,y) - cp e (t,x,y) 
at the point (to, xq, yo), for any e\ > 0, we can find c,d 6 1R and X, y € Sfc, such that: 



(c, *(xo - yo)\x - yo| 27 " 2 + #(2 7 + 2)x \x \ 2 \X) e J 2 > + (u io (t , x )), 

(-d, ^(x - y )[x - yol 27 " 2 - ^(27 + 2) 2/0 |y | 27 ,y) e J 2 -"((l - A)^ (t , yo)), 

c + d = D t <p e (t , x , y ) = -f3 - 4 and finally 



-(i + ||B(to,aro,lto)||)/< 



X 



< B(t ,x , yo) + eiB(t , x , yo) 2 



-Y 

Taking now into account (15. 5D . and the definition of viscosity solution, we get: 

-c- iTr[a*(t ,zo)Xa(t ,zro)] - (^(x - y )|x - y | 27 ~ 2 + 

6{2 7 + 2)xo\x \ 2 ^b{to,x )}-^ lo {t ,xo) < and 

d- lTr[(T*(t ,y )Ya(t ,y )} - (*(x - y )\x - yo| 27 ~ 2 - 

0(27 + 2)y o |yo| 27 ,^o,yo)) - (1 - A)^ (to,yo) > 



(5.8) 
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which implies that: 

-c-d < ±Tr{a*(t ,xo)Xa(t ,x ) - a*(to,yo)Ya(to,y )} 
+ (j(xo - yo)\x - yo\ 2j ~ 2 ,b(t ,x ) - b(t ,y )} 
+<0(2 7 + 2)x |x | 27 , b(t ,x )) + (0(2j + 2)y (2/o| 27 , b(t , y )) 
+tpi(t ,x ) - (1 - A)^j(t ,yo)- 
But from ()5.7[) there exist two constants C and C\ such that: 



(5.9) 



(x ,y )\\ < C\x - yol 27 " 2 and (||a 2 (x )|| V ||a 2 (y )||) < C x 6. 



As 



1/ ai(a?o,yo) -ai(x ,yo)\ a 2 (x ) 
B = B(t ,x ,y ) = - | + 

-ai{x ,y ) ai(x ,y ) I \ a 2 (yo) 



e 



then 

^ < 

£ , . ... \_j j 

It follows that: 



B<-\xo-yo\ 2 ' ( - 2 ( 1 ' i + C X 9L 



5 + eijB 2< C( I| Xo _ yo |2 7 -2 + £l |xo _ yo| 4 7 -4 ) ^ ^| +c , i0/ (51Q) 

where C and Ci which hereafter may change from line to line. Choosing now ei = e, yields the relation 
B + e x B 2 < £(|:c - y | 2 ^ 2 + |x - yo| 47 " 4 ) ^ ^ | + Ctfl. (5.11) 
Now, from (|2TT|) . (IBTHI) and (l57TT|) we get: 

^Tr[(r*(to, x )Xa(t , x ) - <r*(t , 2/o)*V(t , yo)] < ~(\x - y | 27 + |x - yo| 47 " 2 ) + CiO(l + |x | 2 + |y | 2 ). 

— ■ € 



Next 



7 C 2 
(-(»o - yo)\xo - yo\ 2j ~ 2 ,b(t ,x ) - b(t ,y )) < —\x - yo| 27 



and finally, 

(0(2 7 + 2)x |x | 27 , &(*o,*o)> + (^(27 + 2)y |y | 27 , H*o,yo)) < OCQ. + |x | 27+2 + |y | 27+2 )- 

So that by plugging into (|5.9p and note that A > we obtain: 

< ?(ko - yol 27 + \x - y | 47 " 2 ) + Ci0(l + kol 2 + |yo| 2 ) + ?N - yo| 27 + 
ec(i + |x | 27+2 + M 27+2 ) +4>i (to,x ) - (i - A)^ (to,yo). 

By sending e— >0, A— >0, >0 and taking into account of the continuity of ipi and 7 > 2, we obtain 
[3 < which is a contradiction. The proof of Theorem 2] is now complete. □ 

As a by-product we have the following corollary: 
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Corollary 1 Let (v , ...,v m ) be a viscosity solution of {2.J$ which satisfies a polynomial growth condi- 
tion then for i = 1, m and (t, x) £ [0, T] x ]R k , 

rT 

v\t,x)= sup E[ i> Us {s,Xf)ds-^ guT UTn {T n ,X%)\ [Tn<T] ]. 

mm Jt n>l 

Acknowledgement: the authors thank gratefully Prof. J.Zhang for the fructuous discussions during 
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